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ABSTRACT 

Knowing the amount of ionizing photons from young star-forming galaxies is of particular impor- 
tance to understanding the reionization process. Here we report initial results of Subaru/Suprime-Cam 
deep imaging observation of the SSA22 proto-cluster region at z = 3.09, using a special narrow-band 
filter to optimally trace ionizing radiation from galaxies at z ~ 3. The unique wide fleld-of-view of 
Suprime-Cam enabled us to search for ionizing photons from 198 galaxies (73 Lyman break galax- 
ies (LBGs) and 125 Ly-a emitters (LAEs)) with spectroscopically measured redshifts z ~ 3.1. We 
detected ionizing radiation from 7 LBGs, as well as from 10 LAE candidates. Some of the detected 
galaxies show significant spatial offsets of ionizing radiation from non-ionizing UV emission. For some 
LBGs the observed non- ionizing UV to Lyman continuum flux density ratios are smaller than values 
expected from population synthesis models with a standard Salpeter initial mass function (IMF) with 
moderate dust attenuation (which is suggested from the observed UV slopes), even if we assume very 
transparent IGM along the sightlines of these objects. This implies an intrinsically bluer spectral 
energy distribution, e.g, that produced by a top-heavy IMF, for these LBGs. The observed flux desity 
ratios of non- ionizing UV to ionizing radiation of 7 detected LBGs range from 2.4 to 23.8 and the 
median is 6.6. The observed flux density ratios of the detected LAEs are even smaller than LBGs, if 
they are truly at z ~ 3.1. We find that the median value of the flux density ratio for the deteced LBGs 
suggest that their escape fractions is likely to be higher than 4%, if the Lyman continuum escape is 
isotropic. The results imply that some of the LBGs in the proto-cluster at z ~ 3 have the escape 
fraction significantly higher than that of galaxies (in a general field) at z ~ 1 studied previously. 
Subject headings: galaxies: evolution — galaxies: high-redshift — cosmology: observations — inter- 
galactic medium — diffuse radiation 



1. INTRODUCTION 

Ionizing radiation from star-forming galaxies is a likely 
primary source of cosmic reionization. Although the 
ratio of the flux density of Lyman continuum escap- 
ing from a galaxy to that produced in the galaxy, the 

1 Based on data collected with the Subaru Telescope, which is 
operated by the National Astronomical Observatory of Japan. 

2 Some of the data presented in this paper were obtained from 
the Multimission Archive at the Space Telescope Science Institute 
(MAST). STScI is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS5-26555. 

3 Okayama Astrophysical Observatory, National Astronomical 
Observatory of Japan, Honjo, Kamogata, Asakuchi, Okayama 719- 
0232; iwata@oao.nao.ac.jp 

4 College of General Education, Osaka Sangyo University, 3-1- 
1 Nakagaito, Daito, Okaka 574-8530, Japan; akinoue@las.osaka- 
sandai.ac.jp 

5 Department of Astronomy, Graduate School of Science, Kyoto 
University, Kyoto 606-8502, Japan 

6 Optical and Infrared Astronomy Division, National Astronom- 
ical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588; 
yuichi.matsuda@nao.ac.jp 

7 Subaru Telescope, National Astronomical Observatory of 
Japan, 650 North A'ohoku Place, Hilo, HI 96720.; furu- 
sawa@subaru.naoj.org 

8 Research Center for Neutrino Science, Graduate School 
of Science, Tohoku University, Sendai 980-8578, Japan; 
haya@awa.tohoku.ac.jp, kousai@awa.tohoku.ac.jp 

9 Astronomical Institute, Graduate School of Science, To- 
hoku University, Aramaki, Aoba-ku, Sendai 980-8578, Japan; 
akiy ama@ ast r . tohoku .ac.jp, yamada@ ast r . tohoku .ac.jp 

10 Laboratoire d'Astrophysique de Marseille, UMR 6110, 
CNRS/Universite de Provence, 38 rue Joliot-Curie, 13388 
Marseille Cedex 13, France; Denis.Burgarella@oamp.fr, jean- 
michel.deharveng@oamp.fr 



so-called escape fraction (/ eS c)> is a key parameter for 
evaluating the contribution of galaxies to cosmic reion- 
ization, it has been poorly constrained due to the fact 
that Lyman continuum photons are easily absorbed by 
the intergalactic medium (IGM). Direct observations of 
Lyman continuum from galaxies at z > 4 are virtu- 
ally impossible because of a rapid increase in the num- 
ber dens ity of Lyman limi t systems toward high red- 
shifts (|Inoue fc Iw ata 2008). Therefore, we must focus 
on z ~ 3 where the IGM opti c al de pth is still about 
unity on average. ISteidel et al.l ([2001) detected Lyman 
continuum in a composite spectrum of 29 Lyman break 
galaxies (LBGs) at (z) — 3.4. The observed flux density 
ratio (i.e., without correction for IGM absorp tion) was 
/1500//900 = 17.7 ± 3.8. iShaplev et all (|200a hereafter 
S06) detected Lyman continuum individually from two 
out of 14 observed LBGs at z ~ 3 through a deep long- 
slit spectroscopy, and derived /1500//900 = 12. 7 ±1.8 and 
7.5 ± 1.0 for the two LBGs. For an average of 14 objects, 
/1500//900 = 58 ± 35. So far the number of galaxies with 
direct observation of Lyman continuum is too small to 
reveal properties shared by galaxies with large / esc , or to 
estimate a typical value of / esc . 

In this paper we report the initial results of the obser- 
vations of Lyman conti nuum from galaxies a t z ~ 3 with 
Subaru/Suprime-Cam (|Mivazaki et al.ll2002j ). A narrow- 
band filter imaging was a dopted to search for Lyman con- 
tinuum, as pioneered by llnoue et"aTI (|2QQ5h for galaxies 
at z ~ 3. The use of narrow-band imaging instead of 
slit spectroscopy enables us to examine ionizing radiation 
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from a large number of galaxies simultaneously, as well as 
to examine emission offsets from rest-frame non-ionizing 
UV radiation. The target field of the present study is the 
SSA22 field where the prominent pr oto- cluster of galax- 
ies at z = 3.09 has been dis covered ([Steidel et al.lll998L 
12000b lHavashino et al.ll2QQ4l . hereafter H04). We adopt 
a cosmology with the parameters = 0.3, = 0.7 
and Ho = 70 km/s/Mpc. All magnitudes are measured 
in the AB system. 

2. OBSERVATIONS AND DATA REDUCTION 

We have built a special narrow-band filter NB359 with 
a central wavelength of 359 nm and a FWHM of 15 
nm (and > 10% transmittance between 350 nm and 371 
nm). This filter is designed to trace Lyman continuum 
of galaxies at z > 3.06 with a contamination from non- 
ionizing radiation less than 1% for a typical star- forming 
galaxy. Laboratory measurements have shown that the 
transmission of the filter is less than 0.01% in the wave- 
length range of 400 nm-1200 nm. We also confirmed 
that the central wavelength is stable with rms=0.47 nm 
at every position on the filter. In Figure [T] we show the 
transmission curve of the NB359 filter as well as trans- 
mission curves of NB479 which was used in the survey 
for Ly-a emitters (LAEs) at z = 3.1 (H04), V, R, i' 
band filters and spectra of model star-forming galaxies 
at z = 3.09, which is used in section 4.3. 

The imaging observations using Subaru/Suprime-Cam 
with the NB359 filter were carried out in 2007 Septem- 
ber 10-14 (UT). Three nights were photometric while in 
the other two nights cloud coverage was relatively high. 
The seeing condition was good throughout the observ- 
ing dates, with a FWHM of 0."5-0."9 in NB359 images. 
The pointing center (a = 22 h 17 m 26 s .l, 6 = +00 o 16'll" 
[J2000]) is aligned to the previous Subaru/Suprime-Cam 
observations of the SSA22 field (Matsuda et al. 2004, 
hereafter M04; H04). For half of the shots the position 
angle was switched from 90°to 270°, in order to reduce 
the effect of chip-to-chip sensitivity variation. 

The data reduction was made using SD FRED version 
1.2.5 (|Yagi et al.ll2002l : lOuchi et alJl2QQ4 i. The rejection 
of frames under poor sky condition led to a selection of 
45 frames with 1,800 sec exposures to produce the mo- 
saic image. The FWHM of point sources is 0."8. The 
flux calibration was made using measurements of spec- 
troscopic standard stars. The 3 a limiting magnitude 
of the final mosaic image is estimated to be 27.33 AB 
mag for l/'2 diameter aperture. Source detec tion was 
done using SExtractor (|Bertin fc Arnoutslll996l ) version 
2.5.0, using both i^-band and NB359 images. Galac- 
tic extinction has been corrected us ing; the dust map by 
iSchlegel. Finkbeiner fc Davis] (|1998T ). 

3. RESULTS 

There are 198 objects with spectroscopically measured 
redshifts larger than 3.0 in the Suprime-Cam field-of- 
view (32' x 26'). Ten objects classified as QSOs/AGNs 
from optical spectra and/or X-ray observation have al- 
ready been excluded. Among the 19 8 objects there are 
44 LBGs reported in the literature ([Steidel et al.ll2003l . 
S06), 29 LBGs selected as /7-dropout galaxies and iden- 
tified spectroscopically with VLT/VIMOS (Kousai et 
al., in preparation), and 125 Ly-a emitters (LAEs) and 
Ly-a 'blobs' selected through a narrow-band (NB497) 



imaging with Subaru/Suprime-Cam (H04, M04) and 
identified spectr oscopically with Su baru/FOCAS and 
Keck/DEIMOS (jMatsuda et al.ll2006l Matsuda et al., in 
preparation). 

We detected 16 objects with secure detection (> 3<r) 
within l."2 diameter apertures in NB359; 6 are LBGs 
and 10 are LAEs. The redshift range of the LBGs is 
3.04 < z < 3.31, while for the LAEs it is 3.07 < z < 3.10. 
In order to eliminate a possibility of spurious detections, 
we split the frames used to create the final NB359 im- 
age into two sets and generated two mosaic images. All 
16 objects were detected in both images at > 2<j level. 
We also executed a detection test using a negative im- 
age in NB359 and found that a probability of spurious 
detection at > 3a level is ~ 0.4%. Therefore, we are 
confident in the reality of source detection in NB359. In 
Figures [2] and [3] we show 5" x 5" postage stamp images 
of the detected objects. The objects shown in Figure 
[2] are within fields observerd with the Advanced Cam- 
era for Surve ys (ACS) on board the Hubble Space Tele - 
scope (HST; lAbraham et "all l2QQ7t iGeach et all [2007), 
and their F814W images are shown in the figure. 

In addition, the object SSA22a-C49, detected in Ly- 
man continuum by S06, has 2.95<r significance in our im- 
age. The measured flux density within the l."2 aperture 
(5.5 ±1.0 x 10 -2 /iJy) is consistent with that measured 
with a long-slit spectroscopy by S06 (6.9 ± 1.0 x 10 -2 
/xJy). We add this object to the detection list in NB359, 
bringing the total number of detected objects to 17 (7 
LBGs and 10 LAEs). The images of SSA22a-C49 are 
shown in Figure [2j On the other hand, we could not 
find any significant signal in NB359 at the position of 
SSA22a-D3 which is another object detected in Lyman 
continuum by S06. Since the 3<r detection limit within 
l."2 aperture at the position of this object (5.8 x 10 -2 
//Jy) is well below the flux density reported by S06 
(11.8 ± 1.1 x 10 -2 /iJy), it is not clear why the object is 
not visible in our NB359 image. The images of SSA22a- 
D3 are also presented in Figure 03 

All of the 10 LAEs detected in NB359 show a promi- 
nent emission line around 497 nm in their spectra taken 
with FOCAS and DEIMOS, but continua are not well 
detected due to their faintness. All but one spectrum of 
these 10 objects have a resolution high enough to distin- 
guish the doublet if the emission line detected with the 
NB497 filter is [O n] A3727 at z = 0.33, and we find no 
such feature. However, there is still a possibility that 
they are AGNs at lower redshifts if the emission line is 
C IV A1549 at z = 2.21, [C in] A1909 at z = 1.60 or 
Mg II A2798 at z = 0.78, although some of them show a 
spatially extended emission in the NB497 image. Unfor- 
tunately, the wavelength coverage of the spectra obtained 
so far is too narrow to definitely rule out the low redshift 
possibility. Follow-up spectroscopy would be required to 
clarify whether these objects are really LAEs at z c± 3.09. 

If a faint foreground object lies very close to an ob- 
ject at z ~ 3.1, they might mimic Lyman continuum in 
our NB359, and it would be difficult to distinguish it 
at longer wavelengths if the object at z p± 3.1 is brighter 
than the foreground object. ISiana et "all (|2007l ) discussed 
such a possibility of contamination by faint foreground 
objects i n the z y 3 L yman continuum survey. Follow- 
ing |Siana~eFaI] ([2007), we roughly estimate the proba- 
bility of such a case by assuming a surface density of 
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faint galaxies and their spatial distribution to be un- 
corrected with the z ~ 3 objects. The apparent mag- 
nitudes of galaxies detected in NB359 range from 26.5 
mag to 27.5 mag. The surface number density of galax- 
ies in this magnitude ran ge from [/-band number count 
bv IWilliams et al.l (jl996h is - 10 5 mag -1 deg~ 2 . If we 
consider that a foreground object within a l."0 radius 
cannot be distinguished in our NB359 image, each ob- 
ject has ~2.4% chance of such foreground contamination. 
Since we have 198 spectroscopic sample galaxies, several 
objects among the galaxies detected in NB359 may be ex- 
plained by contamination by such a foreground object. 
However, it would be difficult to imagine that all 17 de- 
tections could be due to foreground contamination. 

4. DISCUSSION 

4.1. Spatial Offsets of Lyman Continuum for LB Gs 

As seen in Figures [2] and 03 shapes and positions of the 
emitting regions in the 17 detected galaxies in NB359 
are quite different from those in the i?-band, especially 
for LBGs (marked with circles in Figures [2] and [3]) . The 
rms of positional offsets between the NB359 and the R- 
band images for foreground objects with similar magni- 
tude range (23.5 < R < 27.0) is - 0."25. The offsets 
of 4 LBGs among the 7 LBGs detected in Lyman con- 
tinuum exceed 3<r, and the average offset is 0."97 (3.8<r, 
corresponding to 7.4 kpc at z = 3.09). Higher resolu- 
tion of HST/ACS images clearly shows that in each ob- 
ject with positional offset between NB359 and i?-band 
images there are substructures in F814W which agree 
with the positions of the regions detected in NB359. 
Except for a few possible cases of foreground contam- 
inations (as discussed in section 3), such differences in 
the shape and the position between the emitting regions 
of ionizing radiation and those of non-ionizing UV may 
give us a clue to understand how Lyman continuum es- 
capes from galaxies. For instance, Lyman continuum 
may escape throug h a chimney-like structure in the in - 
terstellar medium (|Razoumov fc Sommer-Larsenll2007h . 
and we may see the emission only from some limited re- 
gions. Another possibility is that the spatial distribution 
of the Lyman continuum sources is different from that 
of the non-ionizing UV emitting stars. As shown in Fig- 
ure El the higher spatial resolution of HST/ACS reveals 
that many of the detected galaxies have multiple knotty 
structures and that the Lyman continuum emission ap- 
pears to be associated with one of such knots. Note that 
at the moment with the F814W image we cannot distin- 
guish between superposition of a foreground object and 
multi ple knotty substructu res frequently seen in LBGs 
(e.g., Lowent hal et al.lll997l ). The position offsets of Ly- 
man continuum from non-ionizing UV peak also imply a 
possibility that long-slit spectroscopy may miss Lyman 
continuum. 

4.2. Photometry and Colors of the Detected Objects 

Because of the various spatial offsets in the Lyman 
continuum, using a fixed aperture size to measure the 
flux for all detected sources is not appropriate to ex- 
amine the flux ratio between the Lyman continuum and 
non- ionizing UV emission from these systems. In order 
to obtain total flux in each image we put the apertures 
with different radii centered at the peak of an object in 



i?-band, examined the curves of growth in NB359 and R- 
band and determined the aperture sizes which are large 
enough to contain the bulk of fluxes in both images and 
which do not contain flux from nearby objects. The 
aperture diameters for the 17 detected galaxies are 2."0- 
4. // 0. In FigureHthe NB359-# colors of the 198 objects 
with spectroscopic redshifts larger than 3.0 are plotted 
against their i?-band magnitudes. The NB359— R color 
corresponds to the apparent flux density ratio of ionizing 
to non- ionizing UV photons for objects at z ~ 3. For 
objects detected in the NB359 image, variable aperture 
sizes described above are used to measure colors, while 
using arrows we also show colors measured with a l."2 
aperture, which was used in object detection, both in 
NB359 and #-band. We find that luminous (R < 25) 
objects, which ro ughly correspond to those w ith L > L* 
of LBGs at z ~ 3 (Sawicki fc Thompsonll2006l ). show rel- 
atively red colors of NB359— R > 2. On the other hand, 
some of the less luminous objects show bluer colors, 
< NB359 -R<2. If the IGM opacity does not depend 
on the source luminosity, this suggests that the ionizing- 
to-non-ioni zing UV escape flux density rat io (a proxy 
of /esc; see llnoue. Iwata fc Deharvend l2QQ6h is lower in 
relatively luminous objects. Interestingly, this appears 
contr ary to the argument by Gne din. Kravtsov. fc Chenl 
(|2QQ8h that / esc is smaller for galaxies with smaller star 
formation rates (or mass). However, we should note that, 
as seen in Figure EJ Lyman continuum appear to be emit- 
ted from substructures in UV-luminous objects, and the 
NB359— R colors of such substructures would be much 
bluer than the colors using total flux densities. 

4.3. Comparison with Model Stellar Populations 

We use deep multi-band optical imaging data (H04) 
to investigate spectral energy distribution (SED) of the 
detected sample. The UV spectral slopes of the detected 
LBGs are flat (0 < V - i f < 0.35). There are two clearly 
distinctive sub-groups of the detected LAEs from their 
V—i' colors: one is red (0.4 < V—i' < 0.8), and the other 
is blue (—0.5 < V — i' < 0). In Figure [H different symbols 
are used for different types of SEDs. Blue LAEs show 
extremely blue NB359-i? colors (< 0.5), while LBGs 
have NB359— R larger than 0.9. The difference between 
these three types becomes clear if we place them into the 
NB359— R and V — i' two-color plane, as shown in Fig- 
ure [5l To compare with these observed colors we calcu- 
lated predictions of colors for young star-forming galax- 
ies with a populatio n synthesis code (Starburst99 version 
5.1: lLeitherer et allll999l ). A model at zero age, with the 
Salpeter IMF (0.1-120M Q ; ISalpeterl fl955h and Padova 
evolutionary tracks with metallicity Z = 4 x 10 -4 is used 
as a fiducial one. Since models with higher metallicity or 
older age have redder colors, this model has the bluest 
SED under an assumption of the Salpeter IMF. The col- 
ors of the model without IGM and dust attenuation and 
assuming / esc = 1 at z = 3.0 and 3.3 are shown as filled 
and open circles connected with a solid line in Figure 
In Figure 03(a) the model does not include flux from neb- 
ular continuum and only flux from stellar sources is con- 
sidered. An arrow indicates the di rection of dust atten- 
uation following a prescription by ICalzetti et al.l ([2000) 
for a galaxy at z = 3. For A < 1200 A we simply extrap- 
olated their attenuation law. Such smooth extrapolation 
toward extreme-UV may be reasonable up to A ^ 800 A 
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(Draine 2003). IGM attenuation has a large dispersion 
for different lines of sight. We show colors with a me- 
dian value of IGM attenuation for z = 3.0 and z = 3.2 
cases, with the optical depth range cont aining 68% of 
all sig htlines of Monte- Carlo simulations (|Inoue fc Iwatal 
2008). The shaded area shows colors of model galaxies 
which can be explained with dust and/or IGM attenua- 
tion. The dashed line in the area indicates colors with 
IGM opacity with 1% probability of occurrence (i.e., very 
rare transparent sight-lines). 

The V — i' colors of the detected LBGs suggest moder- 
ate attenuation by dust, E(B — V) ~ 0.1-0.3, consistent 
with previous studies (e.g., Ilwata. Inoue fc Burgarellal 
2005). However, the NB359-i? colors of the detected 
LBGs are difficult to explain with these models. Even if 
we assume no IGM attenuation, NB359— R colors of the 
three bluest LBGs cannot be explained. One caveat is 
that the model SEDs considered in Figure 03(a) do not 
include nebular continuum. The output SEDs by Star- 
burst99 code have nebular continuum emission assum- 
ing all stellar ionizing photons are converted to nebular 
continuum photons. In Figure [51(b) we show colors of 
the model with both stellar and nebular continua (i.e., 
the model is not self-consistent in regard to photon bud- 
get). The inclusion of nebular continuum makes SEDs 
redder in both NB359— R and V — i' ', and in this case 
the colors of the four redder LBGs would be explained 
with moderate IGM attenuation. However, there are 
still three LBGs whose colors cannot be explained with 
the bluest SED and Calzetti's dust attenuation law. A 
QSO-like spectrum may be able to explain these colors. 
The power-law slope of QSO UV spectrum have large 
variation. Here we assume two-component power law 
(F u oc v a ) with a break at llOOA, as it is made in the 
appendix of llnoue fc Iwat a (2008), and consider the pos- 
sible bluest color of a QSO at z ^ 3. The power law 
index at A < 1100 A can be as steep as —0.6 ([Scott et al.l 
l2004f ). At A > 1100 A t he steepest slope would be 

r 1.8 ([Shang et all 120051 ). In such a case NB359-# 

and V — i' are ~ 0.23 and 0.07, respectively, without 
IGM attenuation. Thus by considering IGM attenuation 
such QSO-like spectrum would be able to explain the 
colors of the LBGs shown in Figure El although these 
objec ts are not classified as AGNs, according to the cat- 
alog (ISteidel et alJl2QQ3h , 

To explain these strong emission in Lyman continuum 
without considering possibility of QSOs, models which 
can produce much bluer intrinsic colors - such as those 
with a top-heavy IMF - would be required. Indeed, 
model SEDs with a top-heavy IMF with a slope 1.0 
(much steeper than the slope of the Salpeter IMF, 2.35) 
have ~ 0.3 mag bluer color in NB359— R. We also ex- 
amined colors of the SEDs with zero metallicity stellar 
population (Population III) wit h the Salpeter IM F (mass 
range from 1M to 100M o ) bv lSchaererl (|2003f ). In Fig- 
ure [5](a) and (b) the colors of the model at z = 3.0-3.3 
without and with nebular continuum are plotted, respec- 
tively. If no nebular continuum is included (Figure 03(a)), 
the NB359-i? colors of the Pop-Ill models are ~ 0.3 mag 
bluer than those of Starburst99 Z = 4 x 10 -4 models 
(with Salpeter IMF) and are close to the top-heavy cases. 

Alternatively, dust attenuation law might be quite dif- 
ferent from that bv ICalzetti et al.l ((2000). For exam- 



ple, if the amount of dust attenuation at 900A is much 
smaller than that at 1500 A, the dust attenuation arrow 
approaches vertical in Figure [5j and the NB359— R and 
V — i' colors of the LBGs may be explained simultane- 
ously. 

Colors of the detected LAEs are much harder to ex- 
plain with the models considered here. If higher mass 
cut-off (5OM -5OOM ) i s adopted , NB3 59-A color of the 
Population III model bv lSchaererl (|2003l ) (using Salpeter 
IMF and no IGM attenuation is considered) is as small 
as —0.5, and it would be compatible with the bluest 
LAEs. Thus we would like to stress that, although some 
of the NB359— R colors of the detected objects are unex- 
pectedly blue, such colors are not physically unexplain- 
able and they can be reproduced if we consider metal- 
poor stellar populations and/or top-heavy IMF. Note 
that there is a suggestion of the existence of metal-free 
stellar populations at z ~ 3 based on observed data 
([Jimenez fc Haima n 2006). For a reference, there are 
three [O n] emitters at z ~ 0.33 found by the follow-up 
spectroscopy of the SSA22 LAE survey (Matsuda et al., 
in preparation). Their colors in NB359— R and V — i' are 
0.9±0.3 and 0.1±0.1, respectively, and they are similar to 
some of LBGs with Lyman continuum detections and are 
different from those of the detected LAEs. Anyway, deep 
spectroscopy with wide wavelength coverage is required 
to confirm their redshifts and explore their nature. In 
the following discussion for constraints on / esc , we only 
use the detected LBG sample. 

4.4. Constraints on the Escape Fraction 

Th e relative escape fraction is defined as (Ste idel et al.l 
l200lh : 



esc,rel 



(/l50o/ /90o)i: 



-exp(n G M. 



900 ) 



(1) 



(/l50o/ /90o)obs 

where (/i50o//90o)int and (/i5oo//90o)obs are the intrinsic 
and observed UV to Lyman continuum flux density ra- 
tios, respectively, and tigm,900 is the IGM optical depth 
for Lyman continuum photons along the line of sight. 
For the ease of comparison with previous studies on Ly- 
man continuum at z ~ 3 (Steidel et al. 2001; Inoue et 
al. 2005; S06), first we use (/i50o//90o)int = 3.0. If the 
dust attenuation at 1500 A A1500 is known, / e sc,rei can be 
converted to / esc as / es c = 10 _0,4Al500 / e sc,rei ([Inoue et al.l 
2005: [Siana et al]l2007l ). We calculate the observed flux 
densities at rest-frame 1500 A by interpolating flux den- 
sities in V and i?-bands, and use flux densities with the 
NB359 filter as / 900 . obs- For different LBGs, different 
aperture sizes are used to eoncompass the bulk of fluxes 
and to avoid contamination from neighbours (see sec- 
tion 4.2). For the seven detected LBGs, (/i5oo//90o)obs 
ranges from 2.4 to 23.8 with a median value of 6.6. If we 
consider the case without IGM attenuation, this median 
value corresponds to f eS c,re\ — 0.45, and it is a lower limit 
of /esc,rei (for (/i50o//90o)int = 3.0). Under the assump- 
tion of dust attenuation E(B — V) = 0.15, which would 
be reasonable to explain their V — i' colors (see Figure |5j), 
a lower limit of / e sc is 0.11. If a correction for the IGM 
attenuation for a median opacity through the NB359 fil- 
ter t igm = 0.59 (for z = 3.0) is applied, / esc = 0.20. 
One unevitable obstacle to estimate escape fraction from 
observed flux density ratios is the uncertainty of the in- 
trinsic UV to Lyman continuum flux density ratio which 
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cannot be observed. For the bluest model galaxy gener- 
ated with Starburst99 (top-heavy IMF, low metallicity, 
no dust attenuation, no nebular continuum) considered 
in section 4.3, (/i5oo//90o)int is 1-07. If we adopt this 
intrinsic flux ratio instead of 3.0, estimated values of the 
escape fraction drop. For the case of no IGM attenua- 
tion and E(B — V) = 0.15, / esc ,rei is 0.16 and / esc is 0.04. 
In Table [1] we summarize our estimates of f esc ,rei and 
/ esc for L BGs, as well a s the a verage values of f esc ,rei ob- 
tained bv ISteidel et all (|200lh and S06. Although there 
are uncertainties on (/i5oo//90o)int? amount of dust at- 
tenuation and IGM opacity, the high f esc ,re\ for the LBGs 
at z ~ 3 we obtained makes a clear contrast with a strin- 
gent upper limits / eS c.rei < - 08 fo r stacked images of 
z ~ 1.3 galaxies (Sia na et al.ll2007l ). Also, five among 
the seven detected LBGs have (/i5oo//90o)obs less than 
10. Thus the UV-to-LC ratios of the LBGs detected by 
us are smaller than th e average values of LBGs obtained 
bv lSteidel et al.l ((2001) and S06. However, we should em- 
phasize that, as it is seen in Figured! these LBGs should 
be the systems with smallest (/i5oo//90o)obs among lu- 
minous galaxies at z ~ 3. The average f esc of the entire 
galaxy populations (at least for L > L* galaxies) should 
be much lower than the value derived for our detected 



LBG sample. The constraints on / esc for UV- magnitude 
limited sample with the present data will be discussed in 
a forthcoming paper. Also, we should remind that the 
SSA22 field is a prominent proto-cluster at z = 3.09, and 
the properties of ionizing radiation in such an uncommon 
environment may not be representative of the epoch. A 
census over a blank field would be required to investigate 
this point. 
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TABLE 1 

Compilation of observed UV to LC flux ratios (/i5oo//90o)obs and escape fractions of LBGs at z ~ 3. 



References 



Sample 



(/l500//900)obs 



(/l500//900)int a 



fei 



Steidel et al. (2001) 
Shau lev et al. (2006) 
Shaplev et al. (2006) 
This work 
This work 
This work 
This work 



29 LBGs, Average 
2 LBGs, Direct 
14 LBGs, Average 
7 LBGs, Direct 
7 LBGs, Direct 
7 LBGs, Direct 
7 LBGs, Direct 



17.7±3.8 
12.7±1.8, 7.5±1.0 
58±25 
6.6 (median) 
6.6 (median) 
6.6 (median) 
6.6 (median) 



3.0 
3.0 
3.0 
3.0 
3.0 
1.07 
1.07 



0.31 
0.43,0.72 
0.094 
0.46 b 
0.83 d 
0.16 b 
0.30 d 



0.11 c 
0.20 c 
0.04 c 
0.07 c 



a (/i50o//90o)int * s an assumed intrinsic flux ratio used to derive / e sc,rel an d /esc from the observed UV to LC flux ratio. (/i5oo//90o)int — 3.0 
has been assumed in the previous studies shown here, and the case (/i50o//90o)int — 1-07 corresponds to the bluest model galaxy SED 
generated with the Starburst 99 code. 
b no IGM attenuation is assumed. 
c E(B -V) = 0.15 is assumed. 

d a median opacity at z = 3.0 through the NB359 filter tigm = 0-59 is assumed. 
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Fig. 1. — Filter transmission curves for the filters used in this study and spectra of a model star-forming galaxy at z = 3.09 generated 
with the Starburst99 code. The solid line represents a model without IGM attenuation, and the dashed line is for a model with average 
IGM attenuation for objects at z = 3. See section 4.3 for details of the model. 
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Fig. 2. — Postage stamp images of objects which are detected in the NB359 image and HST/ACS F814W images are available for them. 
In each panel NB359 (left), i?-band (middle) and F814W (right) images are shown. Field of view is 5" x 5". For the i?-band and F814W 
images a contour map of the NB359 image (2 and 3 a) is over-plotted. Symbols at the lower-left of the NB359 images show the object 
type: filled circles: LBG, filled triangle: 'blue' LAEs. See section 4.3 for the definition of types. The images of SSA22a-C49, which is one 
of the two objects repoted to be detected in S06 and is detected in our NB359 image with 2.95cr level, is labelled. 



NB359:LyCont R:UV(1600A) NB359:LyCont R:UV(1600A) NB359:LyCont R:UV(1600A) 




Fig. 3. — Postage stamp images of objects which are detected in the NB359 image (HST/ACS F814W images are not available). In each 
panel NB359 (left) and i?-band (right) images are shown. Field of view is 5 /; x 5 /; . For the i?-band images a contour map of the NB359 
imag e (2 and 3 a) i s over-plotted. Here we also show images of SSA22a-D3, which is one of the two objects reported to be detected in 
Shaplev et al. ( 2006) but is not detected in our NB359 image. Symbols at the lower-left of the NB359 image show the object type: filled 
circle: LBG, filled triangle: 'blue' LAEs, and filled squares: 'red' LAEs. See section 4.3 for the definition of types. 
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Fig. 4. — NB359— R colors of objects in the SSA22 field with spectroscopic redshifts z > 3.0. The higher and lower dashed lines indicate 
approximately la and 3cr limits in the NB359 image. The objects without significant detection (< 3cr) in the NB359 image are plotted 
with open circles and upper arrows. The detected LBGs are shown with filled circles. LAEs are classified with two types according to their 
rest-frame UV slopes: red UV colors (filled squares) and blue UV colors (triangles). See text for details. For the objects detected in the 
NB359 image, NB359— R colors are measured with variable aperture sizes centered at peak positions in i?-band determined to encompass 
total fluxes, and arrows indicate how colors change if l."2 apertures centered at peaks in NB359 and i?-band are adopted. 
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Fig. 5. — (a) NB359-i? and V - i' colors of objects detected in NB359. LBGs detected in NB359 are shown with circles, and two LAE 
sub-groups divided by their UV slopes are shown with open squares and triangles. Filled and open circles connected with a olid line are the 
color tracks of a model galaxy with the bluest SED with the Salpeter IMF from z = 3.0 and 3.3, respectively. The shaded area indicates 
a color range which can be explained with attenuation by d ust and IGM on this model SED. The arrow indicates the direction of dust 
attenuation following a prescription by Calzetti et al. (2000) and changes in colors with E(B — V) = 0.1 attenuation. Two filled circles 
with horizontal error bars show the colors of the bluest model SED with the median IGM attenuation at z = 3.0 and z = 3.2. The error 
bars represent the ranges of I GM opacity with 68% probability. The expected colors with the stellar population at zero metallicity based 
on the SED by Schaerer (2003) from z = 3.0 to z = 3.3 are also shown (labeled with "Pop-Ill"). The model galaxy SEDs do not include 
nebular continuum emission, (b) same as (a) for observed colors of the objects detected in NB359, but model galaxy SEDs include nebular 
continuum emission. 



